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Abstract 

Utilizing the Higgs Effective Theory, including dimension eight operators, in conjuction with 
^ spinor-helicity and unitarity methods, I present analytic expressions for the virtual corrections for 

gg — > Hg and qq — ► Hg amplitudes at order a^/mf . These contributions become enhanced as the 
invariant final state mass increases. 
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I. INTRODUCTION 

The Higgs boson is the last remaining piece of the standard model to be observed, and 
its observation has been the impetus behind much theoretical and experimental work. The 
dominant production mechanism at the LHC is through a gluon initiated top loop, and this 
mode of production has been well studied to high orders in both the top mass and strong 
coupling expansions. A different but related observable of interest is the Higgs produced 
with a large transverse momentum. This observable has received considerable attention 
[H El E], both a probe of new physics [I] and as experimental avenue for Higgs discovery [3]. 
To produce a Higgs with a large transverse momentum, it is necessary for it to recoil off a 
hard jet. This makes it particularly difficult for study, as the gluon fusion process begins at 
one loop, and higher order QCD corrections to the Higgs pt are two-loop four-point func- 
tions. Practically, one can either tackle the higher order full theory contributions through 
asymptotic expansions P [7J |8], or extend the effective theory, obtained from integrating 
out the top quark, to include more heavy mass effects [9]. Extending the effective theory 
is useful from several perspectives. Corrections in the effective theory are universal, in the 
sense that incorporating new heavy physics is a matter of adjusting matching coefficients. 
Further, by gauge invariance, once one has the matching coefficients, one can include higher 
order heavy mass and loop corrections to many other processes, including the phenomeno- 
logically important Higgs production in association with two jets, relevant for determination 
of Higgs coupling to the weak vector bosons. All one needs is a one loop correction in the 
effective theory, where methods of loop calculations have become quite sophisticated. In this 
paper I present the one loop corrections to Higgs production with a hard jet in the extended 
Higgs Effective theory. With these corrections one can reproduce the full theory amplitude 
to order 0{a 2 /m\), when combined with the matching coefficients for two-loop hard region. 

II. THE BASIS 

The operator basis used in calculating the matching coefficients in an earlier paper [9J , is 



O 




o 2 



HD a F^D a F a ^ 
vM 2 



3 



^a/i Tpbv jpca £abc 
v a fi J 



1 = 



vM 2 



V 



2 



° A ~ ^M 2 ° 5 " ^M 2 (1) 

Where v is the vev of the Higgs particle and M is the heavy scale at which the operators 
are generated. In the Standard Model, this would be the mass of the top quark. 

Using the equations of motion, to this order in the large mass expansion, we can rewrite 
the last two operators in terms of quark bilinears as 

HD a Fg V DpF a P v ^ 2 Hq llfl T a qiq -^T% 

vM 2 ^ 9 vM 2 { ' 

i=X,3=l 

HF^WDPF^ _ ^ , HFg l/ D u q i 'y a T a q i 

vM 2 L> 19 vM 2 [ } 

i=i 

where rif is the number of quark flavors. One can immediately see that the first operator 
will not contribute to any event involving three partons at this order: one needs to have a 
four partonic event for it to contribute, or higher loop corrections. It is further advantageous 
to decompose the operator O2 as 

HD a F* u D a F a ^ u Hd a d a (F^ v F a » u ) HF^ v D u D^F^ a HF^F^F? f abc 
r.\l 2 2vM 2 +4 r.\l 2 1 r.\l 2 
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= -^0 1 + 40 5 -20 3 (4) 

where I have made use of the fact that on-shell the total momentum flowing through the FF 
factor in the first operator must be negative that of the Higgs, if all particles are considered 
incoming. Thus the partial derivatives reduce to an over-all factor of m\. The relation 
itself follows from expanding out one of the partial derivatives in the first operator and then 
making use of the Bianchi identity. With this decomposition, it is necessary to only calculate 
the one loop corrections to O5 and O3, as the Ox amplitudes are well known [TO] . 



III. METHODS 



Typically calculating one loop corrections to higher dimensional operators can be daunt- 
ing due to the high rank tensor integrals encountered. However, with the use of spinor 
helicity methods and generalized unitarity, one can tame the difficulty of such calculations 
in the case of straight QCD and super Yang-Mills theories [UJ H2]. Further, the end results 
and intermediate steps exhibit a simplicity of form generally obscured by traditional methods 
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of evaluation. The simplicity of analytic results come from using gauge invariant and on- 
shell calculational techniques (and the simple basis of master integrals at one loop). As such 
there is no reason to suppose that higher dimensional gauge invariant lagrangians ought to 
yield any more complicated scattering amplitudes. The only subtlies that can arise are from 
the rational parts of the one-loop amplitudes. If one takes strictly four dimensional cuts, 
there are then undetermined rational contributions to the amplitude. Suitable methods for 
determining these rational contributions to amplitudes, such as on-shell recursion, typically 
break down for amplitudes at loop level. Specifically, there is a contribution from infinity 
when taking the large z limit of the amplitude for many choices of the onshell deformation 
of the amplitude. Often this large z contribution can be worked out from auxiliary shifts 
[13] [TH [15]. For higher dimensional operators, this is even more problematic, as they have 
contributions from infinity even at tree-level [16]. In this regard the standard Higgs effective 
theory involving just the HFF operator is the exception, due both its low dimensional order 
and its resemblance to pure QCD [TTl IT81 [HH 120] . 

Thus for my purposes, it is advantageous to make use of d-dimensional unitarity [2H E21 
[23] 124] to calculate the cuts, as by calculating cuts in d-dimensions (that is, the phase space 
of the cut is d-dimensional), one can determine the full analytic structure of the amplitude 
as long as the intermediate states are massless. Also, cutting in d-dimensions, one can still 
make use of spinor integration to perform the tensor reduction [25l |26| [27] . This method is 
more direct than the standard Passarino-Veltmann reduction (particularly for many powers 
of loop momentum in the numerator), as the whole process of reduction to scalar integrals 
amounts to the extraction of residues in the spinor variables. These residues correspond 
directly to the coefficients of scalar integrals. Furthermore, working directly with spinor 
variables, the simple analytic structure becomes manifest. 

IV. CONVENTIONS 

There are only two types of three partonic events in the Higgs production process: those 
involving two quarks (or quark and anti-quark) and a gluon, or three gluon events. For each 
operator we present the helicity partial amplitudes following the conventions of so that a 
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full amplitude has either of the two forms: 




(6) 



(5) 



where C would be the appropriate matching coefficient, and the index i denotes the incoming 
momenta with helicity Aj. Subscripts q or q denote quark or anti-quark momenta, and will 
be suppressed where appropriate. 

The amplitudes are presented here unrenormalized, so that the correct counterterms 
must still be added. I give the minimal helicity amplitudes, other helicity configurations 
either vanish, or are obtainable from parity or charge conjugation. We use Weyl spinors 
|p±) and denote their products as (pj — \pj+) = (ij) and (pi + \pj— ) = [ij] with the 
QCD sign conventions = (i\pj\i] = 2pi ■ Pj = Sij. For spinor chains, one writes 

(i|a|j] = (ia)[aj] = a ll (i\'y IJ, \j] and (i\ab\j) = (ia)[ab](bj) with a and b null vectors. The 
SU (N c ) fundamental matrices are normalized as tr{T a T h ) = 5 ab and [T a ,T b ] = iV2f abc T c . 

I regularize keeping the external states 4-dimensional, and extending the loop momenta 
in 4 — 2e dimensions, and the tensor algebra in 4 — 2eS dimensions. Then one can obtain 
results in t'Hooft-Veltmann scheme j2H] by taking 5 — ► 1 or Four Dimensional Helicity 
scheme [221 ED] with 5^0. 

V. CONSTRUCTING THE CUT INTEGRANDS 

To construct the integrands for the cut amplitudes, one must be careful to make sure that 
they are valid in d-dimensions for the cut momenta. Taking the naive 4-dimensional tree 
amplitudes written in terms of helicity spinors and extrapolating to d-dimensions can lead to 
errors in the rational terms [XT], l3T| 132] . One can then either build tree amplitudes with the 
external states to be used in the loop integration to be d-dimensional, recover the rational 
terms from on-shell recursion or collinear constraints [12] (taking all cuts as strictly four 
dimensional), or build specialized Feynman rules for the rational part [53] • This last course 
of action applied to a non-renormalizable theory requires knowledge of the rational parts 
generated by higher rank tensor integrals, which, while straightforward to obtain, have not 
been classified. Then a matching calculation is needed to extract out the loop level rational 
contributions to the Feynman rules. 



5 



Thus in a theory with a non-renormalizable power counting, one can either use Feynman 
diagrams or an extension of off-shell recursion [31], or try to tackle the on-shell recursion 
with a boundary term. Given the low number of partons, it is easiest to build the integrands 
from Feynman diagrams using color-ordered Feynman rules. Each side of the cut amplitude 
is taken to be the sum of Feynman diagrams in Feynman gauge contributing to the tree 
level process on that side of the cut. Any gluon lines from the tree amplitude running across 
the cut are then truncated, and the polarization sum is replaced by the metric tensor. This 
leaves only the ghost diagrams to be calculated separately. To check gauge invariance of the 
final result, I calculate the cut for several choices of the reference momentum of the external 
polarizations. The final result then is independent of these choices. 

VI. SPINOR INTEGRATION 

One obtains for a general two particle cut: 

At-iooplK-cut = Z / d D L 1 5 + (Ll)d D L 2 5 + (L 2 2 )5 D (L 2 -Li-K) 

AiA 2 J 

AT\L x 1 \ Vl) ... )Vj -L^)At\L- 2 x \ Pk) ... )Vl ,-L- l ^) (7) 

where K — Pi + ... +Pj — — (pk--- +Pi), Aj are the polarizations of the intermediate states, 
and the momenta is considered to be incoming on all amplitudes. Further, one can integrate 
over the L\ momentum in favor of L 2 and write 

At-^K-cut = Z / d D L6 + (L 2 )6+((L - Kf) 

Ai A2 

A l l ee ((L - K) Xl ,...., -L X2 )A% ee {L- x \ (K - L)~ Al ) (8) 

To make use of spinor integration to perform the scalar reduction, one further decomposes 
the integration variable into 4-dimensional (I) and — 2e dimensional, (/i) pieces, L — I + fi. 
(Since the — 2e dimensions are considered space-like, one interprets /x • \i as — /i 2 , with /i 2 an 
explicitly positive quantity. So (L + P) 2 = (1 + P) 2 — /i 2 , where P is any four dimensional 
vector.) As long as the external states are four dimensional (so p, can only form a product 
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with itself) the integration over the angles in the — 2e dimensions may be performed, so that 
the integration measure becomes: 

jflri+/r2\r+//r rs\2\ _ K^Y f .2 / . .2\ -1-e j4T r+ /T 2 ..2\r+/V7 r^\2 ..2<i 



d D L6 + (L 2 )S + ((L - Kf) = i—^- / d/iV) £ <™+(f - /i 2 )<5+((Z - K) 2 - /i 2 ) (9) 

r(-e) J 

As explained in [23J, given a fixed time-like vector K (taken as the cut momentum), one 
can further decompose the 4-dimensional part of the integration as an integration over a 
single parameter z and a null vector I as I = I + zK. 

J d 4 l = J d 4 l5(l 2 )dz(2l ■ K) (10) 

This decomposition is what allows for spinor integration to be used in d-dimensions. The 
integration over the null vector can be written in terms of spinor variables [35], j d 4 l5(l 2 ) = 
J °° tdt J (Idl) [Idl]. Scalar products 21 -P then satisfy 21 -P = t(l\P\l]. The two delta functions 
from the cuts allow one to integrate over the t and z integrals. The \i integral never need 
formally be evaluated, as scalar integrals in the basis have unique analytic forms in fi. 

One is left then only considering the integration in the spinor products j (Idl) [Idl] . Specif- 
ically, a term in the expression for the cut 1-loop amplitude will transform as, 



r(-e)7 J' " '(i\K\i} m -" +2 nUWM 

The Ri and Qj are /i 2 (=d/ ^ru) dependent four momenta derived from Pi and Kj, 

R, = (vT^)P, + "-^-^ (12) 
, Kf - (1 - ^l^n)K j ■ K 

The actual reduction proceeds by using a form of the Schouten Identity to partial fraction 
the spinor products: 

[al] _ \ab] | [ac] 



[bl][lc] [bc][bl] [bc]{lc] 
Thus the spinor products [/•] may be eliminated in favor of (/•) up to a holomorphic anomaly 
(or vice versa) (274 EH] . 

The final spinor integration reduces to evaluating the holomorphic anomaly and the 
residues of the integrand in the spinor variable \l). These residues can neatly be organized 



into the coefficients of the scalar integrals in the decomposition of an amplitude. The 
required formulas for the coefficients can be found in [23] and were derived in 

The scalar integral basis is built from pentagon, box, triangle, and bubble scalar integrals 
defined in higher dimensions. More precisely, one finds the coefficient for an integral of the 
form: 



^ = ^ £ __1__ (15) 

where n < 5, N is some integer, fc, a sum over a subset of the external momentum and 
D = 4 — 2e. One generates these higher dimensional integrals from the spinor integration, as 
the residues from the integration are in general polynomials in fi 2 multiplying the analytic 
expression for the cut integral when N = 0. These powers of fi 2 can then be reinterpreted 
as part of the measure of a scalar integral in another spacetime dimension [211 [37] , giving 
rise to the higher dimensional integrals. 



VII. CALCULATING WITH HFFF 



Here I present the calculation of the all positive helicity gluon one-loop correction to the 
HFFF operator. The contribution from this operator is particularly easy to evaluate as 
the four dimensional cuts completely determine the amplitude, and due symmetry of the 
helicity configuration, a cut in only one channel needs to be evaluated in order to find the 
full contribution. 

The cut in momentum K = pi + p 2 gives the integral for the process as simply: 

I = -J dL 2(L- P2 ) 2 M 2 6 {L)6 « L ~ K ^ ^ 

At this point, the integral is simple enough that one could just evaluate by using the 
Cutkovsky rules to convert it back into a loop integral, and either Passarino-Veltmann 
reduce, or reduce via Feynman parameters. One then evaluates the cut of the resulting 
expressions. 

However, for illustration purposes, I go ahead and use the full machinery of the spinor 
integration. Following the above steps for setting up the spinor integration, one finds 

6 t 12 ] 2 w Wr 1 - /w[^i^L (i|fluW[ ^,^ |ifa,lq[231 d7) 



M 2 r( _ e) y -v-/ y s (l\K\l} 2 (l\Q\l] 
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where Q = -{y/l^u)p 2 - (1 - v / T Tr ^)f , Ry = (Vl^u)^ + (1 - V±^)Pij • and 
p ij = At thi s point one can directly apply the formulas given in [25J. 

A. Triangle Coefficient 

There is only one triangle integral contributing to the amplitude in this cut. The coef- 
ficient itself is given by the residue of the poles of (l\KQ\l). These poles are easily found 
to be given by the Weyl spinors of the momenta P + = — u)p\ and P_ = — — u)p 2 - 
That is \P+) = |1),|P+] = {y/T^u)\l), and |P_) = |2),|P_] = -(y/T=u)\2]. The residue is 



given by 



C Tri = - lim 6 ^ , 1 r j-((P+ - tP JP 13 Q|P + - rP_) [13] 

+ (P+ - rP_|P 23 Q|P+ - rP_)[23] + (P+ <-> P_)) (18) 
= -<12>[12]™M (19) 



B. Bubble Coefficient 

The bubble coefficient is also straightforward to find. Here there is a auxiliary null vector 
i] introduced when evaluating the holomorphic anomaly. Often there can be subtlies when 
T) is fixed such that it coincides with another pole in the integrand. However, here one can 
take it to be p 3 with no issues and proceed normally. 



6[12] 2 2 V -K (l\R 13 K\l)[13] + (l\R 23 K\l)[23] 
"" M 2 [ (v\K\ri] (l\QK\l) 

1 (P + \ V \P + ](P + \([13]R 13 + [23}R 23 )K\P + ) 
JT^u> (P + \K\P + ]{P + \ V K\P + ) 
(P_MP_](P_|([13]P 13 + [23]P 23 )K|P_> 



(P-\K\P-](P-\rjK\P-) 
12[13][23][12] 

: M 2 



)) (20) 

(21) 
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FIG. 1: A contributing diagram with an insertion of the HFFF operator. 




FIG. 2: A contributing diagram with an insertion of the HFDJ operator. 
C. Full Amplitude 

Combining the bubble and the triangle together, and using the symmetric helicity con- 
figuration, one has, 

mW(l + ,2+,3+) = 6[13] ™ 12] (2/ 2 4 ' 2e (pi +P2) + 2 Pl .p 2 l£-*(p 1 ,p 2 , (p 3+Ph )) 

+ 2lt 2e (pi + Ps) + 2pi • p 3 lt 2e (Pii P3, (P2 + Ph)) 

+ 21^^ + P 3 ) + 2^2 • P3/3" 2e (P2, P3, (Pi + Ph))) 

(22) 

where 7,f is the D dimensional n point scalar integral. 

VIII. RESULTS 

For reference I first give the tree level HFF helicity partial amplitudes. 

<-< 1+ ' 2+ < 3+ ) = M%> (23) 

n ,t„(l + ,2 + ,8-) = -pg SI (24) 

'«V(l 5 ",2, + ,3 + ) = l^ (25) 
At tree level we have for the HFFF operator: 



mWlW ) = ^!™ (26) 
mVF(l + ,2 + ,3-) = (27) 
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At one loop I obtain: 

m HFFF (l\ 2+, 3+) = m° HFFF (l+, 2+, 3+)r r iV c ^(( 2 ^ Z ^t )(( Z ^) £ + + ( Z ^) £ )) + 0(e) 

47T e z (l-2e) -5i 2 -5 2 3 -013 

(28) 

mj, fFF (l + ,2 + ,3>0 + O(e) (29) 

The HFDJ operator has no tree level coupling to gluons, and at loop level couples 
through a quark loop, giving 

o. 



m 



1 

HFDJ 



1+ 2 + , 3+) = -m° HFFF (l + , 2 + , 3 + )n /7 ^ + 0(e) (30) 



47T 



m^ DJ (l+, 2 + , 3-) = -m° HFF (l\ 2+, 3") n/ ^(^^) + 0(e) (31) 
where n/ is the number of light quarks. Note that 

^HFF ^ 

the tree level partial amplitude 

for the HFF operator. 

Finally, I give the amplitudes involving two quarks and a gluon. The HFFF has no tree 
level coupling to quarks at this order, and at loop level, 

<FFF^X q X)=ni° HFF (l~ q X q X)N c ^5(-^) + 0(e) (32) 

The amplitudes for the FHDJ operator are 



m HFDJ (l=,2+, 3+) = -<12>[23] 2 (33) 

roWflf ,2+ 3+) = m ° HFDJ (l-,2+,3 + )(^r T )(N c Ui + ±U 2 ) (34) 

_ 1 ( n ( - m h y , c ~ m M ^ 

e(l-2e) l 12 l -5 13 J +l -5 23 J J 

+ 0(6) (35) 

i^Y + ^ + Oie) (36) 



2+,3+) 


— m HFDj(^q J 2q • 


3+) ? r 






— J23 




5*12 5*13 


35+16 




3523 6S23 


18 


u 2 




3 


" e^-S 12 ) 2 


6(1 - 2e) 



Note that 



\2 



rr " ( ml ] fO^le) • (37) 



,47r^ ormNe r(l + e)r(l - e) 

A 

It is straightforward to check that the presented amplitudes preproduce the correct 
anomalous dimensions where known (the HFFF operator) [38] and the correct scheme 
dependence as given in [39J. 
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IX. CONCLUSIONS 



I have presented the virtual partonic amplitudes necessary to calculate the subleading 
terms in the heavy mass expansion at two-loops to the Higgs p t spectrum, by using unitarity 
based calculational methods on operators with a non-renormazible power counting. While 
this is a straightforward generalization of its established use in renormalizable theories, 
the power counting complicates a naive application of the methods. Though such power 
suppressed operators are often negligible for phenomonology, in some sectors they are the 
simplest way of estimating higher loop effects at a higher mass scale. Having machinery 
that allows for simple calculations with such operators is therefore desireable, and unitarity 
and spinor helicity methods provides such methods. In a forth coming publication I will use 
these amplitudes to determine the subleading effects on the p t differential cross-section. 
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